clc,clear

%画三角波形
T = 1;% 周期
w_0 = 1/T;
t = linspace(0, 2*pi, 1000); % 时间向量
x_triangle = sawtooth(2*pi*t/T, 0.5); % 三角波，0.5参数使得波形在-1到1之间

%傅里叶近似
k = 1:2:41;
fourier_x = 0;
for i = 1:1:20
    fourier_x=fourier_x - 8/(k(i)^2 * pi^2) * cos(2*pi*k(i)*t/T);
end
figure(1)
plot(t, x_triangle, 'r'); % 蓝色线表示三角波
hold on
plot(t,fourier_x,'b')